Influence of spin filtering and spin mixing on the subgap structure of I-V 
characteristics in superconducting quantum point contact. 



O 

o 



in 



a ■ 
o ■ 
o , 

u ■ 

^: 

■ 



O 

o 



o 

(N 
(N 

l> 
O 

o 



X 



I. V. BobkovfEl and A. M. Bobkov 

Institute of Solid State Physics, Chernogolovka, Moscow reg., 142432 Russia 
(Dated: February 1, 2008) 

The effect of spin filtering and spin mixing on the dc electric current for voltage biased magnetic 
quantum point contact with superconducting leads is theoretically studied. The I-V characteristics 
are calculated for the whole range of spin filtering and spin mixing parameters. It is found that 
with increasing of spin filtering the subharmonic step structure of the dc electric current, typical 
for low-transparency junction and junction without considerable spin filtering qualitatively changes. 
In the lower voltage region and for small enough spin mixing the peak structure arises. When 
spin mixing increases the peak subgap structure evolves to the step structure. The voltages where 
subharmonic gap features are located are found to be sensitive to the value of spin filtering. The 
positions of peaks and steps are calculated analytically and the evolution of the subgap structure 
from well-known tunnel limit to the large spin filtering case is explained in terms of multiple Andreev 
refiection (MAR) processes. In particular, it is found that for large spin filtering the subgap feature 
at eVk arises from 2fc*'' and {2k ± 1)*'' order MAR processes, while in the tunnel limit the step at 
eVn is known to result from ■n}^ order MAR process. 

PACS numbers: 74.45.-l-c, 74.50.-|-r 



I. INTRODUCTION 

In contrast to the voltage-biased junctions between 
two normal metallic leads, the I-V characteristics of 
superconducting weak links can be highly nonlinear 
in the subgap region eV < 2A, where A is the su- 
perconducting order parameter in the leads. This 
subharmonic gap structure is explained by multiple 
Andreev reflection (MAR) and has been investigated in 
details in various mesoscopic systems including super- 
conductor/normal metal/superconductor (SNS) plane 
j unctions , superconductor / insulator / superconductor 
(SIS) junctions and quantum point contacts 
both theoreticallji^'^^'^^^i^^^iii^i^^ and 
experimentalljJ^ii^iiiii^ii^iSa. The differential con- 
ductance dJ/dV of SNS, SIS and nonmagnetic quantum 
point contacts exhibits a number of peaks at eV = 2 A /n 
in the limiting regimes of short ballistic and diffusive and 
long (in comparison with the superconducting coherence 
length ^) diffusive weak links. In the intermediate regime 
^ ~ d, when the interplay between proximity effect and 
MARs takes place, the conductance subgap structure of 
SNS voltage-biased diffusive junction modifies and has 
an additional maximum at roughly eV ~ A -|- A^. Here 
Ag is a minigap in the equilibrium density of states in 
the normal region arising due to the proximity effeclji^. 
In addition, it was theoretically predicted that in long 
ballistic SNS junctions coherence effects give rise to 
resonant structures in the dc current due to Andreev 
quantization^^. 

On the other hand, ferromagnetic weak links possess 
additional characteristic parameters: spin filtering and 
spin mixing. Spin filtering is the difference between 
transmission probabilities for spin up and down elec- 
trons D| and D^, which, in particular, offers an oppor- 
tunity to generate and control spin currents by driving 



conducting electrons. Of course, this "spin filtering" ef- 
fect takes place in a spin-active junction between two 
normal leads, but the corresponding spin current is lin- 
ear in the bias voltage. Superconducting hybrid struc- 
tures with ferromagnetic weak links give the possibility 
to obtain highly nonlinear I-V characteristics, what can 
be essential for manipulating spin currents in spintronic 
devices. By now there are only several theoretical pa- 
pers devoted to the investigation of voltage-biased super- 
conductor/ferromagnet /superconductor (SFS) quantum 
point contacts. A. Martin- Rodero et.al^ were the first 
to calculate the electric current in SFS quantum point 
contact in the presence of spin filtering effect ^ Di 
and found that the positions of subharmonic gap features 
in the I-V characteristics are shifted as compared with 
the nonmagnetic case and located at eVn ~ A/{\/2n) for 
— 1 and _D| <C 1. Further the influence of another 
generic property of spin-active interface called " spin mix- 
ing" on the voltage biased electric current has been the- 
oretically studied. The essence of "spin mixing" is that 
spin up and spin down electrons acquire different phases 
upon transmission or reflection by a spin-active interface. 
The relative phase shift, so called "spin mixing angle" 
O does not manifests itself in the I-V characteristics of 
the junctions with normal leads, but results in the shift 
of subgap step positions to eVn = A(l -t- cosO)/n in the 
zero-temperature dc electric current of SFS junction even 
without spin filtering^^. Further the simultaneous influ- 
ence of spin filtering and spin mixing on the dc electric 
and spin currents has been analytically studied for the 
case of low-transparency SFS junction D^, <C 1^^. 
It was found that the spin mixing results in the split- 
ting of subgap electric current step positions at non-zero 
temperature leading to eVn = A(l ± cos8)/n. The 
dc spin current also manifests subgap structure of the 
same type in the low-transparency limit, but it shows 
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the odd-even effect: only the steps corresponding to odd 
values of n survive in the I-V characteristics of spin cur- 
rent. It is worth to note that although this analytical 
analysis gives correct positions of the spin current steps, 
the step height is considerably overestimated compared 
to the non-pcrturbative numerical calculation except for 
very tunnel limit D^^i <^ 1^^. 

The most strongly non-linear and rich I-V charac- 
teristics can be obtained in higher transmission SFS 
junctions. Such strongly non-linear I-V characteristics 
for spin current have been numerically calculated very 
recently^^. The present paper is devoted to investigation 
of the dc electric current in the voltage biased SFS quan- 
tum point contact for arbitrary spin filtering and spin 
mixing parameters. For quantum point contact sharp 
subgap features in the electric current only exist when 
at least one of transparencies i'lj is low enough. So it 
makes sense to focus on the two cases: D-^^i ^ 1 and 
-Dj <C 1, ~ 1. We demonstrate how small spin filter- 
ing case continuously evolves to the large spin filtering 
one. It is found that increase of spin filtering qualita- 
tively change the subharmonic gap structure. In particu- 
lar, the step structure of I-V characteristics is changed by 
the peak structure for large spin filtering and small spin 
mixing. It evolves into hump structure with sharp onsets 
upon increasing spin mixing. We analytically study the 
positions of the subharmonic gap features in the limit of 
large spin filtering ^ 1, ~ 1 and their dependence 
on the spin mixing parameter. It is found that the peaks 
and steps are sharp enough for low spin-down transparen- 
cies and under the condition _D| w 1 have the onset volt- 
ages eVn = {ei{Vn) — ej{Vn))/2n, where ei(V) are the 
characteristic energies, corresponding to the poles and 
continuum spectrum edges of the Green's function of the 
system at given voltage. This is in sharp contrast to the 
case of tunnel junction Z)|, <^ 1, where subharmonic 
gap steps take place at eVn = (e-"" — e*"")/n. The dif- 
ference results from the fact that for the tunnel junction 
the n**^ order MAR process opens up new transmission 
channel at voltage V„, while for large spin filtering junc- 
tion the (2n)*'' and {2n+ 1)^^ order MAR processes have 
the onset voltage Vn for spin-up electrons and {2n— l)*-^ 
and {2ny^ order MAR processes for spin-down electrons. 



II. MODEL AND METHOD 

To study the non-equilibrium properties of spin-active 
interfaces in this paper we consider the voltage biased 



one-mode SFS ballistic quantum point contact. It is as- 
sumed that despite the small size of the interface region 
d ^ ^ charging effects can be neglected. The leads are 
supposed to be identical conventional spin-singlet s-wave 
superconductors in the clean limit. 

Our theoretical analysis is based on the non- 
equilibrium quasiclassical theory of superconductivity in 
terms of Riccati amplitudesSS generalized for the case 
of magnetic interfaces22i2^. The fundamental quan- 
tity in non-equilibrium quasiclassical theory of super- 
conductivity is the quasiclassical Green's function g = 
g{pf,R,£,t). It is a 8 X 8 matrix form in the product 
space of Keldysh, particle-hole and spin variables. In 
general, the quasiclassical Green's functions depend on 
space i?, time t variables, the direction of quasiparticle 
Fermi momentum pf and the excitation energy e. In 
our case of one-mode quantum point contact the prob- 
lem is effectively one-dimensional and R = x, where x 
- is the coordinate measured along the normal to the 
junction. The momentum pf has only two values, which 
correspond to incoming and outgoing trajectories. The 
interface is located at cc = 0. 

The electric current should be calculated via Keldysh 
part of the quasiclassical Green's function. For the one- 
mode quantum point contact the corresponding expres- 
sion for the electric current reads as follows 
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T3^o {g'^{pf,x,e,t) - g''^{p^,x,e,t)j 



(1) 



where e is the electron charge. Rq = h/e^ is the quan- 
tum resistance, g^ {pf, x,e,t) is a 4 x 4 Keldysh Green's 
function in the product space of particle-hole and spin 
variables. Pf stands for incoming quasiparticle trajecto- 
ries and p^ for the outgoing ones, fi and ai are Pauli 
matrices in particle-hole and spin spaces, respectively. 

For dealing with interface problems it is convenient to 
express quasiclassical Green's function g in terms of Ric- 
cati coherence functions j^'^ and ^■^^'^^ which measure 
the relative amplitudes for normal-state quasiparticle and 
quasihole excitations, and distribution functions x^ and 
x^ . All these functions are 2x2 matrices in spin space 
and depend on {pf,x,e,t). For definiteness the currents 
are calculated on the left side of the interface. Keldysh 
Green's function for incoming trajectory is parameterized 

by 
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Here subscript 1 means that the corresponding func- 
tions should be taken at a; = —0, arguments pf, 
e, i of all the Riccati functions is omitted for 
brevity. The product (E) of two functions of energy 
and time is defined by the noncommutative convolu- 
tion A (g) S e<(^^^^-^^(^^'>Aie,t)B{e,t). Keldysh 
Green's function g^{p_^) for the outgoing trajectory 
can be obtained from Eqs. ^ hy the sub- 



stitution (7f ,7j;^,xf 



(rf,rf,Xf)(£p and 
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Riccati coherence and distribution functions obey 
Riccati-type transport equations^i^. For clean singlet 
superconductor the equations take the form 
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iidt + Vf^^ — )x^ =7-"® A*(x,t)ia2(X)£ + x^ (g) A(x, t)w2 «> 7 
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Here A(x,t) is the superconducting order parameter. 
The time dependence arises due to non-zero electric po- 
tential which is present at least in one of the supercon- 
ducting leads and can not be removed by the gauge trans- 

'), denoted 



R,A ~R,A ~K 
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formation. The quantities {'j^ ' , 7^^ 
by lower case symbols, are obtained by solving the Riccati 
equations for the appropriate trajectory with the asymp- 
totic conditions, which for spin-singlet s-wave supercon- 
ductor are as follows 
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ifr(e) = (l-|7M^((e-el^,,),i)|2)tanh- 



eVir 



2T 



(6) 



(7) 



where the subscript I, r denotes that the appropriate Ric- 
cati function belongs to the bulk of the left (right) super- 
conductor, while subscript 1(2) stands for functions at 
the left (right) sides of the interface. A is the bulk ab- 
solute value of superconducting order parameter for a 
given temperature, which is assumed to be the same in 
the both superconductors. Vi^r is the electric potential in 
the bulk of left (right) superconductor, so = — V; is 
the voltage bias applied to the junction. Quantities 7/^^'^ 
and if-^ are obtained from Eqs. © and (O, respectively, 
by the operation a(e, t) = a(— e, t)*. 

The superconducting order parameter and electric po- 
tential are assumed to be spatially constant in the su- 
perconductors. Under this assumption the voltage drop 



only occurs at the junction region. These simplifica- 
tions are reasonable for quantum point contact. As 
it can be seen from Eqs. (g]) and ([5]), under the as- 
sumptions above the solutions of Riccati equations for 
{j^'^ (x) , j^'^ (x) , x^ (x) , x^ (x)) do not depend on the 
space variable and coincide with the asymptotic condi- 
tions in the corresponding superconductor. 

The quantities (ff ff , If ), de- 



noted by upper case symbols, are expressed via 

.R,A SR,A ^K^^K^ ^ ^~R,A ~R,A ~K\ 



(7i''"^7r'"^a;i',a; 
elements of the interface scattering matrix S for the 
normal-state electrons and holes with the energies at 
the Fermi surface. The interface iS-matrix is a unitary 
8x8 matrix in the combined spin, particle-hole and 
directional spaces. The explicit structure of iS-matrix in 
directional space is 



S 



Sii S12 
S21 S22 



(8) 



where matrix Su contains spin-dependent reflection am- 
plitudes of normal-state quasiparticles from the inter- 
face in i-th half-space, while Sij with i ^ j incorporates 
spin-dependent transmission amplitudes of normal-state 
quasiparticles from side j. Each element 5*^- is a diagonal 

matrix in particle- hole space Sij — Sij{l + f3)/2 + Sij{l~ 
T3)/2. The most general form of S- matrix for a symmet- 
ric magnetic interface without spin-orbit interaction can 
be written as^^: 



'S'll — ^22 — 



je/2 



-ie/2 



(9) 
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where _R|.| + -D|,x = 1, a = ±1 depending on the par- 
ticular mode!3i. Sij — Sij in the considered problem. 
The particular expressions for {tf',f,r^'.^,X^,X^2) 

terms of {'^^'■^ , x^2t ^^2) 5-matrix elements 
are given in Rc£28. Substituting the Riccati coher- 
ence and distribution functions Eqs. ([6]), ([7]) and 

{f^if,Ti;2^,X^2,X^2) into Eq.®, after some straight- 
forward algebraic manipulations one obtains full Green's 
function for the problem considered. The explicit expres- 
sions are quite cumbersome and the Keldysh part of the 
Green's function is written in the Appendix. The electric 
current is then calculated according to the formula ([T|) via 
the Keldysh part of the Green's function, which is ex- 
pressed by Eqs. (lA.ip and (|A.7|) . When direct voltage is 
applied to the junction, the Green's function and, conse- 
quently, the current are expressed as a sum over harmon- 
ics .g(e,i) = and = Ej-e^™-''*- 

rn m 

We focus on the dc component of the electric current, 
which correspond to m = 0. 

III. RESULTS AND DISCUSSION 

The influence of spin filtering on the I-V characteristics 
of electric current for zero spin mixing angle is demon- 
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These expressions are given for an arbitrary value of spin 
mixing parameter 9. The poles pT|) . of spin-up 

Green's function are identical for left and right sides of 
the interface. The poles of 5o,i(e) £^re obtained from 
Eqs.ini), (ini) by the substitution 6 -6 and adding 
the term eysgna; to the expressions pT|) and p^ . The 
corresponding hole Green's functions at x < have the 
poles at the opposite energies — £1.2^ and — ei,2i- 

With decreasing of D| the poles dTTT) , continuously 
evolve to the position of Andreev bound states at impen- 
etrable magnetic surface £f = sgn [sin(6/2)] cos(0/2) + 
eV{l + sgna;)/32I, and merge with the edge of continuum 
spectrum in the absence of spin mixing = 0. The an- 
alytical expressions for the poles of Green's function at 
-Df < 1 are quite cumbersome, and so we do not write 



strated in FiglTJ The panel (a) shows the evolution of 
subgap structure from low-transparency limit with no 
spin filtering {D^ = Di = 0.1) to large enough spin filter- 
ing when Z)| goes up to unity. It is seen that the subgap 
part of the current-voltage characteristics undergoes con- 
siderable changes. With increasing of spin filtering the 
step structure eV — 2A/ri, typical for SNS structures 
and low-transparency short magnetic junctions modifies 
for eV > A a nd convert s to the peak structure located 
a,t Vn = A/(-\/2n^ -|- 1/2) for eV < A (this expression is 
analytically obtained below). The most pronounced peak 
structure arises for very small values of and smears 
with increasing of as it is seen in FiglTJb). 

We are interested in the case eV ~ A and high enough 
transparencies of the junction for one spin direction. Un- 
der these conditions Keldysh Green's function can not be 
interpreted as a product of density of states and nonequi- 
librium distribution function. There is an external fre- 
quency 2eV in the problem under consideration and the 
above interpretation is only possible in the limit eV <C A. 
Nevertheless the subharmonic gap structure results from 
MAR processes between energies corresponding to the 
poles and continuum spectrum edges of the Green's func- 
tion. Therefore we analyze the characteristic features of 
the Green's function. For the case ^ it can be 
done analytically. The time-independent component of 
electron spin-up Green's function gQ^ij(e) at — > and 
Di = 1 has the poles at energies (see Ea. (IA.16p of Ap- 
pendix) 



(11) 



(12) 

I 

them here. 

To define the subharmonic gap features positions we 
make use of semiconductor representation. Keldysh part 
of electron Green's function g^^ J{—2iTi) for Dj^ ^ is 
presented in FigH] as a function of energy for left and 
right side of the constriction. The temperature is as- 
sumed to have zero value. Let e = correspond to the 
Fermi level of the left superconductor, that is £i?i — 0. 
Then for the right superconductor £f2 = gV . We begin 
by the considering of the tunnel case i'fj — + 0. This 
limit is well-known and we only discuss it for compar- 
ison to the large spin filtering limit. The correspond- 
ing situation is depicted in the panel (a). There are 
no difference between spin-up and spin-down quasipar- 
ticles for 9 = and — > and so spin- up and spin- 
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FIG. 1: (a) Zero temperature I-V characteristics of dc electric 
current for 6 = 0. = 0.1, Dt = 0.1, 0.25, 0.5, 0.75, 1 
from bottom curve to top one. (b) Low-voltage part of I-V 
characteristics for = and _D| = 1. _D| = 0.02 (gray solid 
line), 0.05 (dashed line), 0.1 (dotted line) and 0.2 (solid line). 
The current j is measured in units jRqe/Tc for the panel (a) 
and in units jRQe/lODiTc for the panel (b). The voltage 
\eV\ is measured in units of superconducting order parameter 
A. 



down Green's functions coincide in this case. The only 
characteristic feature of Green's function is a square-root 
singularity at the edge of continuum spectrum. Green's 
functions for left and right side of the junction are de- 
noted by the upper case symbols l,r. In the tunnel limit 
g^{e) — — eV) as it is seen in Fig[2I^a). Let an elec- 
tron ei with energy e — —A pass through the junction 
region to the right. For \eV\ < 2A it undergoes Andreev 
reflection in the right lead and a hole /i2 with energy 
e = A -I- 2eV travels to the left. Then the whole pro- 
cess repeats until an electron in the right lead or a hole 
in the left lead meets continuum spectrum edge. For an 
electron e2k+i the condition is —A — 2eVk = A -I- eV 
and for a hole /i2fe 2eVk + A = —A. Consequently, the 
subharmonic gap steps take place at \eV\ = 2A/n, as it 
should be. 

FiglDJb) represents the large spin-filtering case 
0, — 1. For eV < A the most pronounced characteris- 
tic features of Green's function are pole-like singularities. 
The energies of the poles are given by the expressions 
PT|) and in]). For eV > A the poles merge with the 
continuum spectrum edge and transform to the gap edge 
singularities, which disappear when voltage increases fur- 
ther. The subharmonic gap peaks in this case result from 
MAR processes between the poles of Green's function. It 
is important that the poles of spin-up Green's function 
are located at the same energies for both sides of the 



interface. A spin-up electron ei| with the energy ei| 
passing through the junction region to the right converts 
to Andreev reflected hole with the energy — ei| + 2eV, 
which travels to the left. Then the process repeats until 
an electron e2k+i having the energy ei| — 2eVk meets the 
other pole at 62^, or a hole with the energy — ei| -I- 2eVk 
meets the pole of the hole Green's function at —£21- As 
a result subharmonic gap peaks appear at 



\eVk\ = (e2T -eiT)/2fc 



(13) 



For the absence of spin mixing this equation gives the 
following positions of subgap peaks 



Vk = 



A 



v/2fc2 -I- 1/2 



(14) 



For low enough voltages (in fact, for n > 2) we come 
to the approximate formula by A.Martin-Rodero et.al. 
eVk — A/V2fc21. Analogously, as it is seen from FigO^b), 
for spin-down electron MAR processes repeat until an 
electron e2fc+i having the energy Si^ — 2eVk meets the 
other pole at e2i+2eV, or a hole with the energy — £ix + 
2eVk meets the pole of the hole Green's function at —621- 
It is worth to note that for spin-down Green's function 
the pole energies in the right superconductor are shifted 
to ei_2j. + 2eV in comparison to the pole energies £1,24 in 
the left one. Taking into account the fact that in the left 
superconductor £1,21 — —£2,11 one obtains the same peak 
locations as for spin-up current. However, for spin-up 
current the peak at eVk arises from 2fc"^ and {2k + I)*'' 
order MAR processes, while for spin-down current the 
peak at the same voltage results from {2k— 1)*'' and 2k^^ 
order MAR processes. 

For eV > A the poles of Green's function merge with 
the gap edge and even gap edge singularities disappear 
upon further increase of the voltage. For this reason at 
eV > A MAR processes between gap edges manifest itself 
in I-V characteristic only by the changes of the slope. 

And, finally, intermediate case Z3| ^ 0, D| = 0.5 is 
presented in Fig|2{c). In the intermediate regime Green's 
function manifests two types of singularities: the pole- 
like singularities coming from the large spin filtering case 
and the gap edge singularities, which are characteristic 
for the tunneling limit. When gets smaller the upper 
singularity of tends to A+eV and the pole at £i| disap- 
pears. And, vice versa, for the pole at £21 disappears, 
while the bottom singularity tends to —A. Correspond- 
ingly, spin-up current continuously evolves between the 
two limits. For spin-down current there is an analogous 
picture. 

The subgap structure smears out with increasing of 
Di due to the fact that for non-zero values of Green's 
function poles transform to the broadened peaks. It is 
illustrated in Figllfb), where the subharmonic gap struc- 
ture in the limit = 1 is presented for several values of 
spin-down transparency from = 0.02 up to D| = 0.2. 

Now let us turn to the effect of spin mixing on the sub- 
gap features discussed above. This is illustrated in Fig[3] 
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FIG. 2: Keldysh part of electron Green's function as a function of energy. Energy is measured in units of A. Each figure 
represents Green's functions for left {x — —0) and right {x — +0) sides of the interface. Left column correspond to spin-up 
Green's function and right column to spin-down one. The tunnel limit is presented in panel (a). Panel (b) illustrates the large 
spin filtering case Z)| = 1, and panel (c) corresponds to the intermediate case D| = 0.5, Di — » 0. Spin mixing angle 

= 0. Horizontal arrows show quasiparticles, which take part in MAR processes. Solid arrows describe electrons and dashed 
ones correspond to holes. Hole Green's function is not presented in the figures, but it can be easily obtained from the electron 
one (see text). The particular voltages taken for the figures are: \eV\ = 0.67 for panel (a), \eV\ = 0.34A for panel (b) and 
\eV\ = 0.67A for panel (c). 



for Z?! ^ I, Di = 0.1 and 6 = 0, 0.257r and n/2. As it 
can be seen from Eas. fTTjl . p^ . the voltage range where 
one of the poles exists shrinks with increasing G, while 
it enlarges for the other pole. As a result the peak struc- 
ture, describing by Ea. (|13p only exists at lower voltage 
range when Q grows. For \eV\ > A(l — sin0) one of the 
poles merges with the edge of the continuum spectrum 
and evolves to the gap edge singularity. For this rea- 
son the peaks get wider with increasing spin mixing. At 
the same time the peak onsets are practically not shifted 



due to weak dependence of the pole energies on 0. It 
is seen in FigO for 8 = 0.257r. With further increas- 
ing of spin mixing parameter Q the gap edge singularity 
becomes less pronounced and finally disappears. This 
results in evolving of the peak structure to the hump 
structure with sharp onsets, as it is illustrated in Fig[3] 
for Q — 7r/2. Nevertheless, the hump onsets are well 
described by the equation (fTS]) . 
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FIG. 3: Zero temperature I-V characteristics of dc electric 
current for Df = 1, D| = 0.1. 9 = (solid line), 6 = 0.257r 
(dotted line), Q — n/2 (gray solid line). 



tunnel limit, where the step at eVn is known to result 
from n*'' order MAR process. The spin mixing results 
in shrinking of the voltage range where peak structure 
exists and evolving of the peak structure to the hump 
structure with sharp onsets. 
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IV. CONCLUSION 

The influence of spin filtering and spin mixing, which 
characterize magnetic interface, on the TV characteris- 
tics of dc electric current for magnetic quantum point 
contact with superconducting leads is investigated. It 
is found that with increasing of spin filtering the sub- 
harmonic step structure of the dc electric current, typi- 
cal for low-transparency junction and a junction without 
spin filtering qualitatively changes. In the lower voltage 
region the peak structure arises and for higher voltages 
the slope of I-V characteristic only changes. The subhar- 
monic features are explained in terms of MAR processes 
between Green's function singularities and their positions 
are obtained analytically. In particular, it is found that 
for large spin filtering the subgap feature at eVk arises 
from 2fc"^ and {2k + 1)**^ order MAR processes for spin- 
up current and (2fc— I)*'' and 2k^^ order MAR processes 
for spin-down current. This is in sharp contrast with the 



APPENDIX 

Substituting the Riccati coherence and distribution 

functions Eqs. ©, and {T^^^ ,r^f , X{^^2: ^^2) into 
Eq. ^ , after some straightforward algebra we obtain full 
Green's function. The Keldysh part of the Green's func- 
tion at X — —0 corresponding to the incoming trajectory 
has the following form 

„ = Ki (g) xf (g) - K2 a^i (g Kn 

— 2i7r 

K3 g) if (g — K4 ^"^2^ (8) • (A.l) 

Here * means complex conjugation. ki,2,3,4 are diag- 
onal matrices in spin space and depend on quasiparticle 
energy e and time t. They can be explicitly written in 
terms of S'-matrix elements (O, pUj) as follows 



'^iT.i = 7(e) ^ I (e, t) (g ; 7^77; 1 , 7: (g —r-- 



(A.2) 



^i?^eT*e/2 _ ^2(, „ ^y)^^±^e/2 



(A.4) 



^ v/^eT»e/2 _ ^ ^y^^^^^^^^±^e/2^-2^eVt 
K4t I = ±iK2i-.i <Xi TTT-T^ I , (A. 5) 

I 

where we use a notation Quasiparticle energy e has infinitesimal imaginary part 



/(e-z7A ^)/A , |.|<A 
^^'> ^ 1 (6 - sgne^/^^^)/A, |e| > A . ^^'^^ 
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S > 0. For numerical integration over e we introduce The function A^|(e,<) is inverse to 
finite 6 = 0.003A, which models inelastic scattering. 

Keldysh Green's function ^(Pj) for outgoing trajectory 
can be obtain from Eq. (|A.1[) by the following procedure 



-2iTT 7(e) [ —2iTT 

(7(e))* 



2Re [ki^^] {-V, -e) ^ ] + }> \ . (A.7) 



aoTa(e) = 



i^tae±'®7'(e)7'(e-ey) 



(A.9) 



M a/D^7(e)7(e + ey) 

«lT.i(<^) = ^ , TTT ■ (A. 10) 

/Ri:7eT^0/2 - 72(. + ey)/R;7e±*0/2 



The function is of great importance for us because 
it contains all the essential singularities of the Green's 
function giving rise to the subgap features in the I-V 
characteristics. In general, it has the form 

00 

= E PmTa(e)e''^''"* ■ (A.ll) 

m— — 00 

The condition A~^{e,t) (S) A{e,t) = 1 leads to the follow- 
ing recurrent equation for p„iT,i(e) 

PmTa(e)aoT,i(e + e^m) + 
Pm-iT,i(e - eV)aii^i{e + eV{rn - 1)) + 
/Om+iT,|(e + e^)aiT,i(<^ + eV{m + 1)) = Smo (A. 12) 

with the asymptotic conditions pm{^) ^ for m — > ±cxd. 
In order to obtain Eq. (jA.12p one should use the following 
rule 

ci(e)e2'^^™*®C2(e)e2*^^"* = 
ci(e-eyn)c2(e + eyTO)e2*'=^(™+")* , (A.13) 

which is valid for any integer m and n and can be easily 
deduced from the general definition of the noncommu- 
tative convolution. For arbitrary values of parameters 
D'^ and we find the coefficients Pm'\,i{() numerically 
making use of Eq. (jA.12p . However, for the case ai(e) = 
can be easily found analytically. Then the function 
A does not depend on time and 

A^\{e) = ^ . (A.14) 



There are two special cases when Eq. (|A.14[) is valid. The 
first one is the impenetrable surface — — 0. Then 
the poles of the Green's functions, which are determined 
by the equation 

aoTa(e) = 0, (A.15) 

correspond to the well-known energies of Andreev 
bound states at impenetrable magnetic surface Si i — 
±sgn[sin(e/2)]cos(e/2)27. 

The second limit corresponds to full spin-filtering — 
(or Dij = 0) and arbitrary transparency of the other 
spin channel. In this case the Green's function also has 
the poles. The energies of the poles at a; < can be ob- 
tained from Eg. ljA.lSp . which for spin- up quasiparticles 
is reduced to (we assume to be zero) 

72(e-ey)72(e)e*® = . (A.16) 

The pole equation for spin-down quasiparticles is ob- 
tained from (|A.16[) by the substitution 8 — > — and 
e e + eV. The poles of spin-up Green's function at 
X > are the same as at a; < 0. For spin-down quasi- 
particles they are to be deduced from Eq. (|A.16[) substi- 
tuting e — eV for e and changing sign of Q. The explicit 
expressions of the pole energies are quite cumbersome for 
arbitrary . They are only written for the special case 
£)| = 1 (see Eqs. lfTTI) and (fT2|) ). It is worth to note that in 
the case D^D^ ^ the function A depends on time and, 
consequently, A^^ has no true poles. That is the reason 
for smearing of subgap features when increases. 
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